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In general, the application of nuclear models 
to the calculation of neutron capture and in- 
elastic scattering gamma yields is undertaken 
with the following purposes in mind: 

1. to calculate capture gamma-ray spectra 
where not even partial spectrum meas- 
urements are available. In such cases 
a data vacuum is filled with the under- 
standing that if and when reliable 
measurements appear which seriously con- 
flict with the calculation, the gamma 
yield data will be updated as required. 

2. to fill in gaps in measured capture y- 
ray spectra, generally in the continuum 
region. In such cases available exper- 
imental data are used to aid in the def- 
inition of model parameters. 

3. as a part of data analysis procedures for 
nominally complete capture gamma meas- 
urements such as those done recently at 
GGA for W and U-238. This combining of 
advanced measurement and theoretical tech- 
niques seems to represent the optimum 
approach to the data acquisition problem. 

The gamma-ray cascade model, DUCAL devel- 
oped as part of the Oak Ridge Program to generate 
capture and inelastic scattering gamma yields 
describes the cascade process in terms of para- 
meters which either a) embody statistical assum- 
ptions regarding electric and magnetic multipole 
transition strengths, level densities and spin and 
parity distributions, or b) are fixed by exper- 
iment such as measured energies p spin and parity 
values, and transition probabilities for low lying 
states. Calculation of the latter with an appro- 
priate nuclear model are attempted only for ener- 
gies/mass regions for which either the model has 
been shown to apply or its predictions are judged 
to be closer to the truth than statistical assum- 
ptions. 

With the foregoing limitations in mind, the 
viability of calculated yields may be considered in 


the light of the properties of the nuclei involved, 
experimental information available for the defin- 
ition of model parameters, and the use for which 
the yields are intended. In particular, the DUCAL 
code has been used exclusively to generate gamma 
yields for neutron gamma-ray production cross sec- 
tions to be input to transport codes for the calcul- 
ation of integrated gamma dose at some point in a 

( 2 ) . 

shield. As discussed in an earlier paper, in- 
accuracies in calculated yields which lead to dis- 
crepancies in the calculated differential gamma 
flux do not necessarily result in significant errors 
in the calculated integrated dose. Errors in gamma 
group fluxes which are substantially compensating 
over several neighboring groups tend to be cancel- 
led out in the process of folding the fluxes with 
the flux to the dose conversion ratio. A major 
consideration in the calculation of integrated dose 
is the conservation of gamma energy in the decay 

process. This the cascade model does as a matter 

. ( 1 ) 

of course. 

The fundamental nuclear properties are of par- 
ticular significance where one attempts to calculate 
neutron energy dependent capture gamma yields for an 
element. The number of naturally occurring iso- 
topes, their level structure, and the relative mag- 
nitudes of the binding energies of isotopes formed 
by neutron capture are extremely important. Changes 
in the capture spectrum with neutron energy assoc- 
iated with shifts in the relative magnitudes of cap- 
ture cross sections of isotopes having different 

(3) 

binding energies and level schemes, as in W, can 
often be reasonably accounted for since a) the 
gamma cascade model conserves energy, and b) meas- 
ured low lying level structure is frequently avail- 
able for input to the calculation. In addition, the 
neutron energy dependence of the capture spectrum 
associated with changes in the capture state spin 
and parity resulting in the accessability of diff- 
erent sets of low lying states to El and Ml trans- 
itions can be accounted for to within the approx- 
imation of constant Ml, El and E2 primary trans- 
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Figure 1 Comparison between fitted 
strength calculations of the caputure 
resonance. 
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ition probabilities. An example is seen in the 
calculation of capture gamma spectra for the 35 
keV 3+ and 120 keV 2 resonances in 27 At, shown in 
figures [l] and [2], Differences in the capture 
spectra represented by the dashed curves are 


attributable solely to the nuclear selection rules 
operating in conjunction with fixed Ml, El and E2 
transition strengths. In contrast, changes in 
the capture spectrum with capture state which 
depend upon more subtle wave function character- 
istics clearly cannot be approximated by a model. 

A number of nuclear models have been impl- 
emented to aid in estimating low energy level 
structure and gamma decay schemes used to define 
corresponding parameters in the DUCAL code. These 
include a "unified" model which combines the Nil- 
sson treatment of an odd particle in a deformed, 
rotating potential well with the Faessler-Sheline 




( 4 ) 

excitation , and a model for odd-odd nuclei which 
employs products of single particle Nilsson wave 
functions as basis functions and the Gallagher- 
Moskowski coupling rules for the angular momenta 
Of the odd nucleons. <5) In addition, a parameter- 
ization scheme to extend the applicability of the 

° f t ^ ) foregoing odd ~ A model has been investig- 
ated . This involved an attempt to correlate 
values of a) a spin-dependent inverse moment of 
inertia corresponding to a best fit between meas- 
ured and calculated level energies, and b) a 
deformation parameter to which the E2/M1 ratio is 
sensitive . 
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The foregoing nuclear models have been used 
to generate neutron energy dependent capture gamma- 

^ (3 e 7) S f ° r tungSten ’ tantalum , U-238 and 
iron * . In each case the calculated yields 

have been tested against integral capture gamma 
yield measurements performed at the Oak Ridge 
Tower Shield Facility. 
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